clc;
clear;
close all;

remarks

x=linspace(0,xright,nx);
y=linspace(0,yright,ny);

[xx,yy]=meshgrid(x,y);
[y2,x2]=ndgrid(x,y);


figure(1)

ss=sprintf('pmat.m');
A=load(ss);

subplot(1, 2, 1)
mesh(xx, yy, A');
axis image
axis square
daspect([1 1 0.1]);

ss=sprintf('fmat.m');
B=load(ss);

subplot(1, 2, 2) 
mesh(xx, yy,  B')
axis image
axis square
daspect([1 1 10]);


    




